Shape resonances in nested diffraction gratings 



Angela N. Fantino^, Susana I. Grosz^ and Diana C. Skigin^, 

Grupo de Electromagnetismo Aplicado, 
Departamento de Fisica, 
Facultad de Ciencias Exactas y Naturales, 
Universidad de Buenos Aires, 
Ciudad Universitaria, Pahellon I, 
CI428EHA Buenos Aires, Argentina 
dcs@df.uba.ar 



Abstract 

The diffraction problem of a plane wave impinging on a grating formed by nested cavities is 
solved by means of the modal method, for s and p polarization modes. The cavities are formed 
by perfectly conducting sheets that describe rectangular profiles. The electromagnetic response 
of the grating is analyzed, paying particular attention to the generation of resonances within 
the structure. The dependence of the resonances on the geometrical parameters of the grating 
is studied, and results of far and near field are shown. The results are checked and compared 
with those available in the literature for certain limit cases. 
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1 Introduction 



Current radio detection and ranging (RADAR) technologies employ highly concentrated beams 
of electromagnetic energy to scan an area. The return echo is then processed to detect all 
targets of interest ^J. To prevent this kind of detection, a target can either deflect the beam in 
a direction away from the observer, absorb the energy of the incoming beam or alter the return 
echo by changing its frequency or shape. Giving an aircraft textured skin that could interact 
with an incoming beam adaptively would allow most shapes to be metamorphic to radar. Recent 
research based on sub-wavelength, narrow-band, resonant, corrugated, high impedance surfaces 
have shown promising results in this area However, the speed, bandwidth and adaptability 

of these surfaces are limited. To avoid the narrow-band resonant behavior, lumped structures 
can be replaced by broadband antenna arrays that add a controlled phase delay to the reflected 
signal allowing beam steering away from the observer Fractal planar antennas |S1-|Z| exhibit 
multi-band behavior and are small in electrical size, but their scattering behavior is not widely 
understood. Nested multilayer corrugated surfaces constitute an interesting alternative for this 
purpose. If successful, unwanted radiation due to finite ground planes can be minimized and 
the inter elemental coupling can be reduced, thus obtaining the desired effect. The study of the 
electromagnetic response of the two-layer periodic structure proposed in this paper is expected 
to provide a first approach to the behavior of multilayer surfaces and their advantages as part 
of broadband antennas. 

The resonant behavior of infinitely periodic gratings of rectangular profile was studied by many 
authors, in particular for s polarization [HJ-^Sl- The excitation of surface shape resonances in 
structures comprising cavities have been lately studied for particular profiles of the corrugation. 
The resonant features of an isolated cavity or groove of univalued and multivalued geometries 
have been investigated by means of different implementations of two basic approaches: the 
integral and the modal methods The results show that a strong intensification inside 

the cavity is found for certain wavelengths when its profile is bivalued, such as a slotted cylinder 
or a bottle-shaped groove ^Hl-|2ni- The effect produced by the surface shape resonances on the 
response of an array of cavities has also been studied j21j-|23j. 

The modal method employed here to solve the diffraction problem from an infinite grating 
comprising rectangular grooves and bottle-shaped cavities, was first formulated to deal with 
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simple geometries such as rectangular IP, triangular or semicircular j^S]- Later it was 
generalized for arbitrary shapes of the grooves PHI HZji and recently it was also applied to 
bottle-shaped cavities PJ. The implementation of the modal method for perfectly conducting 
structures results in a simple and efficient way of calculating the external and the internal fields, 
without the need of any sophisticated algorithm. 

The purpose of this paper is to analyze the electromagnetic response of a two-layer perfectly- 
conducting grating comprising grooves and cavities. Particular attention is paid to the resonant 
behavior of this structure and its effect in the field distribution inside and far from the corruga- 
tions. The study of such gratings constitutes a first approach to the analysis of N-layer fractal 
gratings, which are expected to be useful for the control of the phase delay in the scattering 
response of broadband microwave antennas. 

In Sec. 2 we pose the problem and give the details of the modal method applied to the present 
structure, for both fundamental polarization modes. In Sec. 3 we give numerical results for the 
reflected efficiency as well as for the near field. We analyze the dependence of the resonances 
on the geometrical parameters of the structure, and compare the results obtained in limit cases 
with those found in the literature. Finally, some conclusions are summarized in Sec. 4. 

2 Theoretical formulation 

We consider the diffraction of a plane wave by a two-layer nested grating. The structure is 
one-dimensional and infinitely periodic, and each period is formed by a bottle-shaped cavity 
and rectangular grooves, bounded by perfectly conducting sheets (see Fig. 1). 
The wavelength of the incoming wave is A and the angle of incidence is ^o- Since the structure 
and the fields are invariant under traslations in the z direction, the problem can be separated into 
two scalar problems corresponding to the basic modes of polarization: s (electric field parallel to 
the rulings) and p (magnetic field parallel to the rulings). In the upper layer, the period has two 
rectangular grooves of width an and height hi (zones 11 and 13, see Fig. 1), and a central neck 
of the same depth but of width ai2 (zone 12). The second layer is a rectangular cavity which 
occupies the whole period, i.e., its width is d and its depth is /12 (zone 21). The surrounding 
medium is vacuum. In what follows, both polarizations will be treated simoultaneously, denoting 
by / the z-component of the electric field in the s-case and that of the magnetic field in the 
p-case. 
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Figure 1: Sketch of the nested diffraction grating. 

In an homogeneous medium, the rectangular geometry of the structure ahows a separable solu- 
tion of Helmholtz equation in cartesian coordinates. The modal method consists in expanding 
the fields inside the corrugations in their own eigenfunctions (modes) that satisfy by themselves 
the boundary conditions at the sides of the cavities. In the homogeneous region {y > 0), the 
total field / is given by: 

oo 

/(x,y) =e*(°o^-^o^)+ i2^e*(°"^+^"2') , q = s,p (1) 

n=— oo 

where qq = ksinOo, Pq = kcosOo, 

an = ksinOn = ao + n— , (2) 

a 



Pn 



— if k^ > al 



(3) 



i ^Ja\-k'^ if k^ < 

k = 2'k/\ is the wave number and Rn is the unknown complex amplitude of the n — th diffracted 
order . The superscript q denotes the polarization case. 

As mentioned above, the fields inside the corrugations are expressed in terms of infinite series. 
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Each zone ij {ij = 11, 12, 13 or 21) has its own expansion as follows: 

oo 



where 



m=0 



sm 



cos 



an 



for (7 = s 
for q = p 



«12m(^) = ' 



sm 



cos 



0.12 



{x - an) 



-{x - an) 



sm 



cos 



for g = s 

for q = p 
for g = s 



^(x - (an + 012)) forg=p 
for q = s 
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ai 



sm 



cos 
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for q = p 

^fim sin + /ii)] for s polarization 

^iim cos [;Uiim(y + /ii)] for p polarization 



W?2m(y) = [^12m sin (^12™^) + S?2m COS (/ii2my)] 



9 = S,^ , 



Wl3m(y) 



W21m(y) 



Cf3^ sin [/iiim(y + hi)] for s polarization 
Cf^j^ cos [;Uiim(y + /ii)] for p polarization 
Clim sin [/U2im(y + /ii + /12)] for s polarization 
Cfim cos [/U2im(y + /ii + /12)] for polarization 



A;2 



if A;2 > 



A:2 if < 



(4) 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 



and Cj?-^, ^i2m ^i2m unknown complex amplitudes. Notice that ufjQ{x) = 0, and 
then the sum in eq. Q) starts from m = 1 in the s-case. The functions ulimix) satisfy the 
appropriate boundary conditions at x = and at x = an, the functions Ui2m{x) satisfy the 



boundary conditions at x = an and at x = an + ai2) the functions u'l^^{x) satisfy the boundary 
conditions at x = an + ai2 and at x = d, and the functions ^t2im(^) satisfy the boundary 
conditions at x = and at x = d. On the other hand, the functions ^limiu) ^iid ^i3m(y) satisfy 
the boundary condition at y = —hi, and the function v/2imiy) satisfy the boundary condition 
at y = —{hi + /i2), according to the case of polarization. 

To solve the problem, the fields in each zone are matched imposing the boundary conditions at 
the horizontal interfaces y = —hi and y = 0. Then, expression @ with (jH)) and H12|) is matched 
at y = —hi with the field in the central neck of the first layer (zone 12), given by (0]) with © 
and ()1U|). For < x < an and an + ai2 < x < d, a null tangential electric field is required. At 
y = 0, the fields in zones 11, 12 and 13 are matched with those in y > 0, given by eq. All 
these conditions generate four x-dependent equations, that are projected in appropriate bases 
(either the modal functions or the plane waves) to give an infinite system of linear equations for 
the unknown amplitudes, for each polarization case. The explicit expressions of these equations 
can be found in the Appendix. To find the numerical solution to this problem, we truncate 
the modal series in (jlj and the plane wave expansions in (^, and get a matrix system which is 
solved by a standard numerical technique of inversion. 

3 Numerical results 

In this section we analyze the diffraction by nested gratings for s and p polarization, paying 
particular attention to the surface shape resonances that characterize their electromagnetic 
response. All the results presented have been checked to satisfy energy conservation within an 
error less than 10~^^. The results for limit cases such as a comb grating (/12 0) 28 and 
a bottle-shaped grating (hi 0) 21 have also been verified. Even though the method was 
developed for an arbitrary angle of incidence, in all the results presented we consider normal 
incidence. We analyze the grating response as a function of the depth of the bottom cavities 
h2/d and of the width of the neck 012/ d (Figs. 2-6) and show the dependence of the behaviour 
of these structures as a function of the incident wavelength (Figs. 7-8). 

In Fig. 2 we plot the specular efficiency versus h2/d for s polarization, where: hi/d = 0.1, 
X/d = 0.65 and 012/ d = 0.2 (solid) and 0.4 (dashed). It can be noticed that there are minima 
in the efficiency curves for certain values of /i2. These minima are sharper and deeper for the 
narrower neck of the cavities, and can be associated with surface shape resonances. 
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Figure 2: Specular efficiency vs. h2/d for s polarization. The grating parameters are hi/d = 0.1, 
\/d = 0.65, 00 = 0°, ai2/d = 0.2 (solid) and 012/^ = 0.4 (dashed) 

In Table 1 we compare the resonant depths of the cavities of the nested grating for X/d = 0.65 
(/i2) with those corresponding to eigenmodes of a closed rectangular waveguide of side d, which 
are given by 

h n 

— = — , ,m,n€Z>o . 14 

d ^(2d/A)2 - m2 ' ' ^ ^ 

Notice that the first value of m in (|14j) is 1 for s- polarization (the x- dependent part of the 

electric field modes is a pure sine function, and m = would imply null field), whereas it is for 

p-polarization (the x- dependent part of the magnetic field modes is a pure cosine function, and 

here m = implies uniform distribution in the x direction), as it is explicitly shown in eq. (jSj). 

Besides, due to the symmetry imposed by the normally incident plane wave, only odd values of 

m are allowed for s-polarization, and only even values of m for p-polarization. The number of 

decimal places kept in the tables was determined by the resolution necessary to define correctly 

each minimum in the efficiency curve. The sharper dips require more decimals than the smooth 

ones. 

The resonant depths of the nested grating are close to those of the closed waveguide obtained by 
eq. ()14[). even though the present structure is open. The relationship between the width of the 
dip, the quality of the resonance and the interior field can be better understood by inspection 
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3 
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0.9915 
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1.300 


1.29466 
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1.374 


1.35445 
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1.462 


1.45025 
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1.625 


1.339 
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1.710 


1.6342 
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1.718 


1.70125 


s 





6 


1.950 


1.752 


p 



Table 1: Resonant depths h/d for a perfectly conducting rectangular waveguide and resonant 
depths h2/d found in the nested structure for au/d = 0.2, hi/d = 0.1, for a wavelength \/d = 
0.65. 

of the near field. 

In Figs. 3 and 4 we plot the electric field magnitude (relative to the incident field) for different 
resonant depths of the structure considered in Fig. 2. In all the contour plots presented (Figs. 
3, 4, 6 and 8) the black represents maximum intensity, and the gray scale is maintained in all 
three figures. Figs. 3a to 3f correspond to the resonant configurations for ai2/d = 0.2, whereas 
Figs. 4a and 4b correspond to the 1** and 5*^* resonant cases for a\2/d = 0.4. It can be observed 
in Fig. 3, that the contour plots of the inner field in the resonant configurations are similar 
to the configurations expected for the eigenmodes of a rectangular waveguide. Each one of the 
figures can be associated with a certain mn mode of the closed waveguide: Fig. 3a. corresponds 
to the 11 mode. Fig. 3b to the 12, Fig. 3c to the 13, Fig. 3d to the 14, Fig. 3c to the 31, 
and Fig. 3f to the 15. Besides, the interior field is strongly intensified: the ratio between the 
maximum value inside and outside the structure varies between 6 (Fig. 3a) and 35 (Fig. 3e) for 
resonant configurations, whereas out of the resonance the field inside and outside the structure 
has roughly the same value. Notice that the vertical scale is not maintained in all the figures, 
and therefore the depth of the cavities look equal although each one corresponds to a different 
depth. When the neck of the cavities is widened, the quality of the resonances becomes lower. 
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Figure 3: Electric field intensity distribution for a grating with hi/d = 0.1, ayild = 0.2, Xjd = 
0.65, 6*0 = 0° and s polarization (solid line in Fig. 2): (a) h2/d = 0.32675; (b) h2/d = 0.66985; 
(c) h2/d = 1.01305; (d) /la/d = 1.35445; (e) /12/d = 1.45025; (f) h2/d = 1.70125. 

as observed in the dashed curve of Fig. 2. At the same time, the resonant depths shift to lower 
values, thus moving further from the predicted depths for the closed waveguide. This behaviour 

is expected since this structure is less similar to the closed one. 

The contour plots of electric field corresponding to this situation are shown in Fig. 4. Fig. 4a 
corresponds to the first dip {h2/d = 0.23) and Fig. 4b corresponds to the 5th. dip {h2/d = 
1.4335). In the first case there is no intensification: the field takes the same values inside and 
outside the structure. The second case corresponds to the narrower dip of the dashed curve in 
Fig. 2, and consequently to the better resonance for the range of /i2 considered. In this case there 
is an intensification, but it is significantly smaller than that of Fig. 3e (for the same resonant 
mode but a^/d = 0.2 instead of 0.4, and /i2 very close to that of the prevoius case). When the 
minimum of specular efficiency is better localized, the enhancement of the field increases and so 
does the quality of the resonance. 
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(a) (b) 

Figure 4: Electric field intensity distribution for a grating with hi/d = 0.1, ai2/d = 0.4, X/d = 
0.65, 00 = 0° and s polarization (dashed line in Fig. 2): (a) h2/d = 0.23; (b) h2/d = 1.4335. 

The equivalent curve to the solid one in Fig. 2 but for p polarization is shown in Fig. 5. It 
can be noticed that the minima are not as deep and sharp as in the s-case, and this fact is 
also reflected in the near field plots of Fig. 6, where we show the magnetic field for the first 
(/i2/d = 0.037) and the 7th. {h2/d = 1.29466) minima of Fig. 5. 

Not all the minima for the p-case are located close to those of the closed waveguide, listed in 
Table 1. This suggests that in general, these minima are not associated with resonances of 
the lower cavities in the nested grating, as it can be seen, for instance, in Fig. 6a, where we 
observe that the interior field is not intensified. However, there are certain depths that seem to 
fulfill the conditions for a resonant mode. This is the case of the 7th dip, whose corresponding 
magnetic field plot is shown in Fig. 6b. According to Table 1, this dip is located very close to 
that of the waveguide for the mode 04, and it can be clearly observed in the contour plot that 
the field distribution is that corresponding to this mode: nearly uniform distribution in the x 
direction and four half wavelengths in the y direction. Even though this is the most intense 
case, the intensification ratio is of about 6, roughly the same value of the less intensified mode 
of s polarization. 

We have also analyzed the grating response as a function of the neck width ai2, and wc have 
found that for each configuration there is an optimum width for which the specular efficiency 
has its minimum (not shown). This value arises from a compromise relationship between the 
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Figure 5: Specular efficiency vs. h2/d for p polarization. The grating parameters are hi/d = 0.1, 
ai2/d = 0.2, X/d = 0.65, 6*0 = 0°. 

narrowness of the neck and the necessary size of the aperture that allows the field to get into 
the cavity. 

In Figs. 7-8 we explore the response of the nested structure as a function of the wavelength. The 
specular efficiency for both polarization modes is shown in Fig. 7, for a grating with hi/d = 0.1, 
/12/d = 0.32675, ai2/d = 0.2 and = 0°. For X/d greater than one, there is only one diffracted 
order, and the specular efficiency is equal to unity (not shown). For smaller wavelengths, certain 
dips start to appear, that become more frequent as A ^ 0. It is interesting to notice that there 
are certain wavelengths in which the efficiency is nearly zero for one polarization, and at the 
same time it is maximum for the other. For instance, at X/d = 0.65 almost all the specularly 
reflected wave is p-polarized, whereas at X/d = 0.5, the reflected wave is s-polarized. Then, 
this structure behaves as a polarizer for certain wavelengths, which can be varied by properly 
designing the cavities. 

Using eq. H14|) to find the resonant wavelengths of a given waveguide, we get the values listed 
in Table 2, where these values are compared with those corresponding to the dips of the nested 
grating with hi/d = 0.1 and ai2/d = 0.2. In the range of wavelengths considered, there are 
five dips for s polarization and eight for p polarization. As it can be noticed in the efficiency 
curve, some of them are sharper than others. This suggests that their corresponding resonances 
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Figure 6: Magnetic field intensity distribution for the grating and incidence conditions of Fig. 
5: (a) h2/d = 0.037; (b) /12/d = 1.29466. 

are also better, in what concerns to the interior field intensification. This can be confirmed by 
inspection of the interior field, shown as contour plots in Fig. 8, for p polarization. The figures 
are ordered by increasing resonant wavelength, and it can be observed that in most of the 
cases the magnetic field distribution is associated with a resonant mode mn of the rectangular 
waveguide. For instance, Fig. 8a is the 02 mode, Fig. 8b is the 60 mode, Fig. 8c is the 40 mode, 
Fig. 8f is the 21 mode, etc. It is evident that the presence of an aperture in the cavity modifies 
the field distribution of each mode, but most of them can still be identified. An interesting 
feature found in the wavelength dependence of the response, is the possibility of having modes 
with n = 0. These modes are not allowed in s polarization since they imply null interior 
field; however, in the p-case, these modes represent an uniform distribution of the field in the y 
direction. This is the case of the modes 20, 40 and 60, shown in Figs. 8h, 8e and 8b, respectively. 
In general, the higher intensification is found when the resonant wavelengths are closer to those 
of the waveguide. Another interesting phenomenon is the splitting of the 41 mode: there are 
two neighbour dips close to the resonant wavelength predicted for the waveguide, but neither of 
them is too close (see Table 2). Even though the interior field configurations corresponding to 
those wavelengths (Figs. 8c and 8d) are different, both can be associated with the 41 mode, but 
none exhibits a significant intensification. The same effect appears in s polarization for the 51 
mode (the corresponding interior field plots are not shown). 
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Figure 7: Specular efficiency vs. \/d for a grating witli hi/d = 0.1, ay^jd = 0.2, h2/d = 0.32675, 
^0 = 0°, for s (solid) and p (dashed) polarization. 

In general, the s-surface shape resonances in partially closed cavities are more significant than 
the p-ones, and this is in agreement with previous reports on multivalued structures [1^00] . 

4 Conclusion 

The diffraction problem of a nested grating with rectangular cavities was solved for both polar- 
ization modes by means of the modal approach. The response of the grating has been analyzed, 
paying particular attention to the surface shape resonances. It was found that these resonances 
are stronger for s than for p polarization: the drops in the specular efficiency are sharper and 
deeper, and the electromagnetic field inside the cavities is significantly intensified. The resonant 
depths and the contour plots of the near field corresponding to the first modes of the structure 
are associated with the eigenmodes of a rectangular waveguide. 

This study is a first approach to the analysis of the behaviour of nested structures in the presence 
of electromagnetic waves. The advantages of using this kind of structures as parts of broadband 
antennas is now being analyzed, not only with the perfectly periodic model but also considering 
a finite structure and multilayers. 
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Figure 8: Magnetic field intensity distribution for the grating of Fig. 7, for p polarization: (a) 
\/d = 0.32099; (b) \/d = 0.331; (c) X/d = 0.35844; (d) \/d = 0.44101; (e) \/d = 0.507; (f) 
X/d = 0.55895; (g) X/d = 0.63671; (h) X/d = 0.88201. 

5 Appendix 

By application of the boundary conditions at the horizontal interfaces, we get four x-dependent 
equations that after appropiate projections generate an infinite system of linear equations for the 
unknown amplitudes. After a little manipulation of the equations, and making the corresponding 
truncations, we can summarize the system for each polarization mode as follows. 

Wa:^ = w''w , (15) 
TO^=r+^V, (16) 
fV = w*+fV , (17) 
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0.55895 
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0.500 


0.507 
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0.466 


0.47762 
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0.397 


0.35844-0.44101 
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0.341 


0.33709-0.357 
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6 





0.333 


0.331 


p 
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0.326 


0.32099 


p 


3 


2 


0.293 


0.30156 


s 



Table 2: Resonant wavelengths for a perfectly conducting rectangular waveguide Xu,/d and for 
the nested grating \g/d, for h2/d = 0.32675, a^/d = 0.2 and hi/d = 0.1. 



where 
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+ 7m) + W4(l - 7i2jfc) for g = s 
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- 
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012021 



/q" cos (^a;) e*"*^(ia; for p polarization 

^ Em JkmJnm^^2\m \tly^^Z iOT q = S 
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Jq"^^ sin (y^x^ sin (^^(a; + aii)^ dx for s polarization 
Jq"^^ cos (^f^a;^ cos (|^(a; + aii)^ dx for p polarization 
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(25) 



(26) 



(27) 



(28) 



(29) 



In the above expressions the unknown modal amplitudes in the neck of the cavities have been 
redefined: 
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2i^l2m^ ror — 6 

L 2i ^i2m^ ror g — p 

2t 12m^ lOr g — A 

Mi2"i f>1 p-iM2m/»l for /7 — Tl 

2i -"l2m^ lor q — p 



(30) 



(31) 



and TZ'i are the unknown Rayleigh coefficients. 
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